The Crooks relation in optical spectra - 
universality in work distributions for weak local quenches 
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We show that work distributions and non-equilibrium work fluctuation theorems can be measured 
in optical spectra for a wide class of quantum systems. We consider systems where the absorption 
or emission of a photon corresponds to the sudden switch on or off of a local perturbation. For the 
particular case of a weak local perturbation, the Crooks relation establishes a universal relation in 
absorption as well as in emission spectra. Due to a direct relation between the spectra and work 
distribution functions this is equivalent to universal relations in work distributions for weak local 
quenches. As two concrete examples we treat the X-ray edge problem and the Kondo exciton. 
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Equilibrium thermodynamics provides the framework 
for the description of the equilibrium properties of macro- 
scopically large systems. This includes the properties 
of systems in equilibrium states given fixed external pa- 
rameters. Additionally, equilibrium thermodynamics al- 
lows for the description of transitions between different 
equilibrium states even if the system is not in equilib- 
rium in the meantime. The near-equilibrium proper- 
ties of systems driven only slightly out of equilibrium 
are determined by its equilibrium fluctuations accord- 
ing to linear response theory. Starting in 1997 with a 
seminal contribution from Jarzynski p] , the field of non- 
equilibrium work fluctuation theorems opened up. These 
relate a measurable non-equilibrium quantity, the work 
performed, to equilibrium free energies even if the system 
is driven arbitrarily far away from equilibrium. Most im- 
portantly, it is not required to wait for the system to 
equilibrate. 

Suppose a system is initially prepared in an equilibrium 
state at a temperature T = l/kg/3. If the Hamiltonian 
H(t) of the system changes during a time interval from 
to tf according to an arbitrary protocol, energy is not 
conserved and work is performed on the system. In order 
to determine the work done on a system two energy mea- 
surements are necessary leading to the notion that work 
is not an observable [5] . The work ui rather is a random 
variable with a probability distribution function [5] 

P P {u) = J ^ e iut G(t), G(t) = (emo)t e -m a (t f )t^ 

(1) 

Here (• ■ •) denotes the thermal average over the 
initial equilibrium configuration and Hu{tf) = 
U'(tf)H(tf)U(tf) is the time-evolved final Hamiltonian 
with U(tf) the time-evolution operator obeying the dif- 
ferential equation id t U(t) — H(t)U(t). 

Let Pb (w) denote the probability distribution function 
for the backward protocol. Then the Crooks relation, 
first shown for classical systems [5] and later extended to 



closed as well as open quantum systems 4 
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establishes a universal connection between the forward 
and backward processes that only depends on the equi- 
librium free energy difference AF of the final and initial 
state independent of the details of the protocol. The 
Jarzynski relation [JJ is contained in Eq. ^ as a special 
case. 

Experimental tests of the Crooks relation and of its 
variants have been performed in recent years. Among 
these are folding-unfolding experiments of small RNA- 
hairpins where the free energy difference between the 
folded and unfolded state has been extracted using the 
Crooks relation 0[S]. Moreover, the Crooks relation has 
been verified in electrical circuits [9], mechanical oscilla- 
tors [10], for small colloidal particles [UJ and for non- 
thermal systems, specifically a driven two-level system in 
a non-equilibrium steady state |12) . 

In this work it will be shown that work distributions 
and therefore the non-equilibrium work fluctuation the- 
orems such as the Crooks relation can be measured in 
a large class of quantum systems. Namely, optical spec- 
tra of systems where the absorption and emission of a 
photon corresponds to the sudden switch on or off of a 
local perturbation are in fact work distributions. Among 
these are X-ray spectra of simple bulk metals as well as 
of quantum dots coupled to fermionic reservoirs. In such 
systems, the absorption spectra A(u) and emission spec- 
tra E(ui) are proportional to forward and backward work 
distributions such that they are connected to each other 
via the Crooks relation. For the particular case of a weak 
local perturbation, they satisfy the universal relations 



A(w + AF) 
A(-w + AF) 
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(3) 



in second order renormalized perturbation theory. Here 
AF is the free energy difference between the system with 
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and without local perturbation. Notice that an inde- 
pendent measurement of AF is not required to establish 
Eq. ([3]) in an experiment. Actually, Eq. ^ permits a 
determination of AF by plotting the spectrum, see be- 
low. Due to the correspondence between spectra and 
work distributions, Eq. Q implies universal relations for 
work distributions of weak local quenches: 



P f {uj + AF) 
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Here, Pp{ui) is the work distribution for a protocol where 
the local perturbation is suddenly switched on and Pb (oj) 
is the work distribution for the backward protocol. 

Consider a system that is coupled to a monochromatic 
light field of frequency uj. This coupling is assumed to be 
small such that it can be treated perturbatively. In the 
following, we will concentrate on the class of systems in 
which the absorption or emission of a photon corresponds 
to the sudden switch on or off of a local perturbation. In 
the remainder, H denotes the Hamiltonian with the local 
perturbation and Hq without, respectively. Such systems 
have already been discussed extensively in the literature. 
In the X-ray spectra of simple metals, for example, a 
system of free fermions has to adapt to a suddenly created 
or annihilated local potential scatterer [13H15) . Quantum 
dots coupled to fermionic reservoirs can be initialized in 
such a way that the absorption of a photon implements 
the sudden switch on of a Kondo impurity [16] . 

Crooks relation in absorption and emission spectra. 
First, we concentrate on the absorption case, the related 
emission spectra will be discussed below. The absorption 
spectrum for incident light of frequency u in second order 
of the system-light field coupling (Fermi's golden rule) is 
related to a dynamical correlation function via Fourier 
transformation 



A(u) = k a I ^ e*-* G A (t). 



(5) 



The constant ka contains parameters depending on the 
experimental details such as the intensity of the incident 
light beam or the system-light field coupling. Its precise 
value, however, is irrelevant for the derivation of Eq. ([3]). 
We consider those systems where the dynamical correla- 
tion function Ga(£) appearing in Eq. ^ is of the struc- 
ture 



G A {t) = 



1 

~Z~A 



Tr(e-^ Ha e iHot e- 



iHt\ 



Za = Tr (e 
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This is the case for X-ray edge type problems such as the 
original X-ray edge problem or the Kondo exciton. Thus, 
GA(t) is the characteristic function of a work distribution 
for a quench from Hq to H, cf. Eq. ([!]). This is the origin 
of the connection between optical spectra and the work 
distribution that allows to observe the Crooks relation 



in an optics experiment. Recently, X-ray edge singular- 
ities have been observed in work distributions for local 
quenches in an Ising chain at criticality |17j . 

The emission spectrum E(w) corresponding to the 
same setup formally resembles the absorption case, only 
the roles of Hq and H are interchanged and the sign of 
w is altered. The rate of photon emission E(ui) is given 
by: 



E{w) = k e 



dt 
2^ 



*G B (t) 
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with 



G E (t) = -Llr (e-P H e lHt e~ lHot ) , Z E = Tr (e"^) . 
Ze 

(8) 

Hence, E(—uj) is proportional to the work distribution 
for a protocol where the local perturbation is switched 
off, that is precisely the backward process to absorption. 
A direct application of the Crooks relation in Eq. ^ 
therefore yields 



E(lj) 



K A 
( 

k e 



-AF) 
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as an exact result. This relation depends on experimental 
details through the parameters ka and ke- The linear 
scaling of \u(A(uj)/E(uj)) with the frequency ui of the 
light beam, however, is universal with a slope /3. 

Two different measurements are necessary to explore 
this relation in experiment, one to determine the absorp- 
tion and one to obtain the emission spectrum. However, 
the Crooks relation can also be measured in a single ex- 
periment for weak local perturbations: the emission and 
absorption spectra are directly related to each other such 
that in second order renormalized perturbation theory 
Eq. Q holds. This has the additional advantage as op- 
posed to the exact relation in Eq. ^ that also the ex- 
periment specific constants ka and ke drop out. 

Crooks relation in a single spectrum. Suppose V is 
the unitary transformation that diagonalizes the Hamil- 
tonian H. As it is well known for generic weak cou- 
pling impurity problems [18] the diagonalized Hamilto- 
nian, VHV^ — Hq + AF, is given by the unperturbed 
one in the thermodynamic limit up to a constant AF |19j . 
The shift AF is the free energy difference between H and 
Hq. More details about the unitary transformation will 
be given below. This relation between the diagonalized 
Hamiltonian H and Ho allows to write the dynamical cor- 
relation functions GA/Eif) m terms of overlaps between 
V and 0: 



G A {t) 
G E (t) 
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For all the relevant cases, it is possible to represent the 
unitary transformation V as an ordered exponential 



V = Oexpfx] 



(11) 



where \ is antihermitian, x ~ ~X: an( i O denotes 
some ordering prescription. Thus, the adjoint is given 
by V* = Oexp[— x] where O is the ordering prescription 
opposite to O. In the X-ray edge problem no ordering 
prescription is needed and x = i(l + <?)<^(0) with g the 
strength of the potential scatterer and <fi{x) a bosonic 
field |15j . For generic weak coupling problems such as the 
Kondo model at nonzero temperature analyzed later, the 
flow equation approach [IS] provides a general prescrip- 
tion for the construction of the unitary transformation V 
as an ordered exponential of its generator rj(B) 



V = T B exp 



dB r)(B) 



(12) 



where the generator is determined by a set of differential 
equations. Here, the symbol T# orders an r)(B) left of an 
Tj(B') if B > B' as the usual time ordering operator. In 
view of Eq. (flO| and Eq. (11 1, the dynamical correlation 



functions G A /E{t) can be written as the expectation val- 
ues of a single exponential with an appropriate ordering 
prescription Pa/e- Here, only the absorption case will 
be shown in detail, the emission case can be obtained 
analogously, 



G A (t) = (Pa exp [ X - X {t)])e 



iAFt 



(13) 



where (• • •) denotes a thermal average with respect to Hq. 
The ordering prescription Pa orders all time-dependent 
x's to the left of all time-independent ones, the time- 
independent x's according to O, and the time-dependent 
ones according to O. For expectation values of ordered 
exponentials such as in the equation above, there exists 
a general identity [30] 



G A {t) = exp (P A exp [ x - *(t)] - 1) 



-iAFt 



(14) 



relating them to the exponential of a cumulant average 
indicated by an index c. At this point we perform a 
cumulant expansion up to second order in x whose va- 
lidity depends on the strength of the local perturbation 
as we will discuss below. The first cumulant vanishes as 
X can be chosen normal ordered relative to the initial 
state. For the X-ray edge problem the cumulant expan- 
sion up to second order is exact within the validity of the 
bosonization technique as x is linear in bosonic opera- 
tors. For more complicated problems such as the Kondo 
exciton the diagonalizing unitary transformation can be 



obtained by the flow equation framework, see Eq. (12 1 



In this case, the generator rj(B) and thus the operator x 
is proportional to the strength of the local perturbation 
such that in the case of a weak local perturbation the ex- 
pansion is controlled by a small parameter. For systems 



where renormalization effects are important, not only the 
bare coupling has to be small but also the renormalized 
one. 

Implementing the same procedure for the emission case 
G_g(i) one observes that G A (t) and Gs{t) are directly 
related to each other within the cumulant expansion up 
to second order via G A {t)e lAFt = G E {t)e~ lAFt . For the 
spectra this result implies: 



E{uj + AF) 



K A 



A(-oj + AF). 



(15) 



Plugging this relation into the Crooks relation, see 
Eq. (p}]), one directly proves the main result, Eq. pi), in 
second order renormalized perturbation theory. 

In the remainder of this letter, we will discuss two ex- 
amples for the Crooks relation in absorption spectra: the 
X-ray edge problem and the Kondo exciton. 

The X-ray edge problem. In the X-ray edge problem 
the absorption of a photon is accompanied by the sudden 
creation of a local potential scatterer in a sea of noninter- 
acting fermions |14) . Hence, we have H = *^2 k ctc^: 

and H = H(g) = H + (2Tr/L)gJ2kk ,:c k Ck ' : - We consider 
a linearized dispersion Sk — Vpk and set vp = 1. The 
Fourier transform of the absorption spectrum is given 
by H2| 



1 

Za 



S(t) = ^TT[e-P Ho e Hot iP(0)e 



^(0)) (16) 



that is yet not in the desired form as in Eq. ©. Us- 
ing the bosonization technique, the fermionic fields tp(x) 
can be represented in terms of bosonic ones, 4>(x), via 
ij)[x) = or 1 / 2 Fe~ l ^ x ^ with a -1 an ultra-violet cut- 
off [21]. The Klein factor F commutes with H(g) and 
does not contribute to S(t) due to its property FF^ = 1. 
The bosonization identity allows to regard the fermionic 
fields as a unitary transformation acting on H(g) such 
that S(t) cx GA(t)e~ lAt with a constant energy shift A 
that can be absorbed into a redefinition of the threshold 
frequency and Ga(0 is in the desired form: 



GJt) = _Ltt ( e -PH 0e i Ho t e -mi+9)t) . 

Za V / 



(17) 



The diagonalizing transformation V of H(l + g) equals 
V = e i ( 1 +9)^(°) [15] . Although the effective strength of 
the scatterer 1 + g is not small, the cumulant expansion 
stops at second order as the operator in the exponent is 
linear in bosonic operators. Hence, in the range of va- 
lidity of the bosonization treatment, the Crooks relation 
in Eq. ^ holds exactly for the X-ray edge absorption 
spectrum. Comparing bosonization [15 with the exact 
treatment [M], it yields the correct result up to second 
order in g. This restriction stems from the linearization 
of the free fermionic spectrum [15] , 

In view of Eq. ( 16 ) it is an interesting question if cor- 



relation functions in other systems can be transformed 
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Figure 1: Absorption spectrum of a Kondo exciton for differ- 
ent temperatures as a function of the light frequency uj + AF. 
As a reference, a NRG curve for T = 100 Tk is shown taken 
from Tiireci et al. |16| . In the regime \u>\ > T where the NRG 
data is accurate the agreement is excellent. The inset exem- 
plifies the validity of the Crooks relation in the absorption 
spectrum of the Kondo exciton. 



we refer to Ref. [IB]. For frequencies \uj\ > T, however, 
where the NRG data is accurate we observe excellent 
agreement with the results of the flow equation formal- 
ism. The inset shows the validity of Eq. ([3]). The ratio 
A(oj + AF)/A(-lj + AF) is the universal function 
independent of any details. 

Conclusions. We have shown that work distributions 
and thus the non-equilibrium work fluctuation theorems 
can be measured in optical spectra of quantum systems 
such as the X-ray edge problem or the Kondo exciton. 
For weak local perturbations, the Crooks relation estab- 
lishes a universal relation within a single spectrum, ab- 
sorption or emission, cf. Eq. (pj). Due to the correspon- 
dence between optical spectra and work distributions this 
translates directly into universal relations for work dis- 
tributions of weak local quenches, see Eq. Q. 
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into the form as in Eq. (|6j. This may allow to deter- 
mine similar relations or to gain an additional point of 
view by realizing the equivalence to a work distribution. 
One example is the local density of states in a Luttinger 
liquid. 

The Kondo exciton. Recently, Tiireci et al. |16j pro- 
posed an experimental setup for a quantum dot where 
the absorption of a photon corresponds to the sud- 
den switch on of a Kondo impurity. Hence, we have 
H o = J2ka £ k'- C la c kcr- and H = H + J2kk' Jkk'-S ■ s kk >:. 
For details about the Kondo problem see for example 
Ref. [22]. The dynamical correlation function Gvi(t) for 
the absorption spectrum is given by Eq. ^. The di- 
agonalizing unitary transformation V can be obtained 
by the flow equation approach 18J, cf. Eq. (12 1, with 
V(B) = J2kk' ( £ k - £fc') Jkk'(B): S ■ s kk >: in 1-loop order. 
The couplings Jk k /(B) are determined by a set of dif- 
ferential equations |23| . Importantly, the flow equation 
framework includes all the renormalization effects such 
as the emergence of a low-energy scale Tk, the Kondo 
temperature. The absorption spectrum is obtained via 
the cumulant expansion up to second order in the cou- 
pling strength. Its validity is restricted to weak coupling 
problems such that we have to require T ^> T K [2~4"j . 
A plot of the absorption spectrum is shown in Fig. [T] 
for different temperatures. As a reference, a NRG-curve 
for T = 100 Tk obtained by Tiireci et al. [16] for an 
Anderson impurity model in the Kondo regime is in- 
cluded in this figure |25j . In the vicinity of the main 
peak at small \w\ < T, the NRG calculation contains 
an unphysical double peak structure. For more details 
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